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Abstract. Most recent MaxSAT algorithms rely on a succession of calls 
to a SAT solver in order to find an optimal solution. In particular, several 
algorithms take advantage of the ability of SAT solvers to identify unsat- 
ishable subformulas. Usually, these MaxSAT algorithms perform better 
when small unsatisfiable subformulas are found early. However, this is 
not the case in many problem instances, since the whole formula is given 
to the SAT solver in each call. 

In this paper, we propose to partition the MaxSAT formula using a 
resolution-based graph representation. Partitions are then iteratively joi¬ 
ned by using a proximity measure extracted from the graph representa¬ 
tion of the formula. The algorithm ends when only one partition remains 
and the optimal solution is found. Experimental results show that this 
new approach further enhances a state of the art MaxSAT solver to op¬ 
timally solve a larger set of industrial problem instances. 


1 Introduction 

Many real-world problems in different areas such as fault localization in C pro¬ 
grams, design debugging, upgradability of software systems, among other, can be 
solved using Maximum Satisfiability (MaxSAT) [2,9,11,14,24]. The increase us¬ 
age of MaxSAT for problem solving results from the improvements of MaxSAT 
technology in recent years. In the last decade, several new techniques and al¬ 
gorithms have been proposed that improved on previous MaxSAT solvers by 
several orders of magnitude. Moreover, the developments in the underlying SAT 
technology, namely identification of unsatisfiable subformulas and incrementality 
have also been a factor in the improvements of MaxSAT solving. 

MaxSAT solvers for industrial instances are usually based on iterative calls 
to a SAT solver. Moreover, most of these MaxSAT algorithms take advantage of 
the ability of SAT solvers to identify unsatisfiable subformulas. However, in most 
cases, algorithms deal with the whole formula at each call of the SAT solver. As 
a result, unnecessarily large unsatisfiable subformulas can be returned at each 
SAT call, resulting in a slow down of the MaxSAT algorithm. In this work, we 


try to avoid this behavior by partitioning the formula and taking advantage of 
structural information obtained from a formula’s graph representation. 

In this paper, we improve on the current state of the art MaxSAT solving 
by proposing a new unsatisfiability-based algorithm for MaxSAT. The new al¬ 
gorithm integrates several new features, namely: (1) usage of resolution-based 
graphs to represent the MaxSAT formula, (2) partition of soft clauses in the 
MaxSAT formula using the referred representation, (3) usage of structural infor¬ 
mation obtained from the graph representation to drive the merge of partitions 
and, (4) integration of these features into a new fully incremental algorithm 
that improves on the best non-portfolio solver from the last MaxSAT Solver 
Evaluation on several industrial benchmark sets. 

The paper is organized as follows. Section 2 formally defines MaxSAT and 
briefly reviews the MaxSAT algorithms more closely related to the proposed ap¬ 
proach. In section 3, graph representations of CNF formulas are described. More¬ 
over, the adaptation of resolution-based graphs is proposed. The new MaxSAT 
algorithm is proposed in section 4. Besides a detailed description, we show how 
to extract structural information from the graph representations and integrate 
it in the new algorithm. Section 5 presents the experimental results of the new 
MaxSAT solver on a large set of industrial benchmark sets used at MaxSAT 
evaluations. Finally, the paper concludes in section 6. 


2 Preliminaries 


A propositional formula in Conjunctive Normal Form (CNF), using n Boolean 
variables Xi, X 2 , ■ ■ ■ ,Xn, is defined as a conjunction of clauses, where a clause is a 
disjunction of literals. A literal is either a variable Xi or its complement Xi. The 
Propositional Satisfiability (SAT) problem consists of deciding whether there 
exists a truth assignment to the variables such that the formula is satisfied. 

The Maximum Satisfiability (MaxSAT) can be seen as an optimization ver¬ 
sion of the SAT problem. In MaxSAT, the objective is to And an assignment 
to the variables of a CNF formula that minimizes the number of unsatisfied 
clauses. Notice that minimizing the number of unsatisfied clauses is equivalent 
to maximizing the number of satisfied clauses. 

In a partial MaxSAT formula ip = iphU (ps, some clauses are considered as 
hard {(ph), while others are declared as soft (ips)- The goal in partial MaxSAT 
is to And an assignment to the formula variables such that all hard clauses in 
iph are satisfied, while minimizing the number of unsatisfied soft clauses in ip^. 
There are also weighted variants of MaxSAT where soft clauses are associated 
with weights greater than or equal to 1. In this case, the objective is to satisfy 
all hard clauses and minimize the total weight of unsatisfied soft clauses. In this 
paper, we focus solely on partial MaxSAT, but the proposed approach can be 
generalized to its weighted variants. Furthermore, in all algorithms we assume 
that the set of hard clauses (ph is satisflable. Otherwise, the MaxSAT formula 
does not have a solution. This can easily be checked through a SAT call on (ph- 
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Algorithm 1: Linear Search Unsat-Sat Algorithm 
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Input: ifi = Vs 

Output: satisfying assignment to v 

{vw, Vr, X) ■<— {iph, 0, 0) 

foreach uji £ ips do 

Vh <—Vr U {ri} // ri is a new relaxation variable 

wr a;* U {ri} 

U {tiJR} 

while true do 

(st, R, vc) ^ SAT(<pw U {cnf(E,,6v„ n < A)}) 

if St = SAT then 

1^ returns // satisfying assignment to v 

A ^ A+ 1 


The most recent state of the art MaxSAT solvers are based on iterative calls 
to a SAT solver. One of the most classic approaches is the linear Sat-Unsat 
algorithm that performs a linear search on the number of unsatisfied clauses. 
In this case, a new relaxation variable is initially added to each soft clause and 
the resulting formula is given to a SAT solver. Whenever a solution is found, a 
new cardinality constraint on the number of relaxation variables is added, such 
that solutions where a higher or equal number of relaxation variables assigned 
the value 1 are excluded. The cardinality constraint is encoded into a set of 
propositional clauses, which are added to the working formula [3,12,16]. The 
algorithm stops when the SAT call is unsatisfiable. As a result, the last solution 
found is an optimal solution of the MaxSAT formula. 

A converse approach is the linear search Unsat-Sat presented in Algorithm 1. 
Here, a lower bound A is maintained between iterations of the algorithm. Initially, 
A is assigned value 0. In each iteration, while the working formula given to the 
SAT solver (line 7) is unsatisfiable, A is incremented (line 10). Otherwise, an 
optimal solution to the MaxSAT formula has been found (line 9). 

Observe that a SAT solver call on a CNF formula tpw returns a triple {st, 
V, vc), where st denotes the status of the solver: satisfiable (SAT) or unsat¬ 
isfiable (UNSAT). If vw is satisfiable, then v stores the model found for ipw- 
Otherwise, vc contains an unsatisfiable subformula that explains a reason for 
the unsatisfiability of vw- 

Several of the most effective algorithms for MaxSAT take advantage of the 
current SAT solvers being able to produce certificates of unsatisfiability. Since 
the SAT solver is able to identify unsatisfiable subformulas, several MaxSAT 
algorithms use it to delay the relaxation of soft clauses. An example is the 
MSU3 algorithm [15] presented in Algorithm 2. Observe that this algorithm also 
performs an Unsat-Sat linear search, but soft clauses are only relaxed when they 
appear in an unsatisfiable subformula. 

Although more sophisticated MaxSAT algorithms exist [19], an implementa¬ 
tion of MSU3 algorithm on the Open-WHO framework was the best performing 
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Algorithm 2: MSU3 Algorithm 


Input: ip = i/Ph U (pa 

Output: satisfying assignment to ip 

1 {pw, Vr, A) ■<- {p, 0, 0) 

2 while true do 

3 (st, R, pc) ^ SkT{pw u 

4 if St = SAT then 

5 return r 

6 foreach uji £ (pc H pa) do 

7 Vb t- Vr U {ri} 

8 ojR ■(— uJiU {ri} 

9 pw {pw \ U {cur} 

10 A ■«— A + 1 


Ti < A)}) 

// satisfying assignment to p 

// ri is a new variable 
// Ui was not previously relaxed 


non-portfolio algorithm at the MaxSAT Solver Evaluation in 2014^. One of the 
crucial features for its success relies on the fact that only one SAT solver instance 
needs to be created [16]. Therefore, a proper implementation of MSU3 should 
take advantage of incrementality in SAT solver technology. In this paper, the 
MSU3 algorithm is further improved with structural information of the problem 
instance to solve. 

3 Graph Representations 

In order to extract structural properties of CNF formulas, different graph-based 
models have been previously proposed. For instance, graph representations have 
been used to characterize industrial SAT instances [1] and to improve on the 
performance of MaxSAT algorithms [18]. In this section, we briefly review the 
Clause-Variable Incidence Graph (CVIG) and adapt the use of Resolution-based 
Graphs (RES) [26] to model relations in CNF formulas. Although other models 
exist [1,25,18], in the context of our algorithm for MaxSAT solving, these were 
found to be the best suited. 

In the CVIG model, a weighted undirected graph G is built such that a 
vertex is added for each variable Xj and for each clause oji occurring in the CNF 
formula p. Moreover, for each variable Xj occurring in clause uJi (either as literal 
Xj or Xj), an edge {uJi,Xj) is added to graph G. The edge weight w{uii,Xj) is 
defined as: 

w{uJ^,Xj) = —^ (1) 

|Wj| 

where jw^j denotes the number of literals in clause Wi and I{xj) is defined as the 
incidence function of Xj in soft clauses as: 

/te) = l+ E o (2) 

Xj^OJ A LJ^iPs 

^ Results available at http://www.maxsat.udl.cat/ 
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Fig. 1. Example of Graph Models 


As described in section 2, several MaxSAT solvers rely on the identification of 
unsatisfiable subformulas. In order to capture sets of clauses more closely related 
that would result in an unsatisfiable subformula, we propose to adapt Resolution 
Graphs (RES) to MaxSAT. 

In the RES model, we have one vertex in graph G for each clause € ip. Let 
Ui and uij denote two clauses such that Xk G Wi and Xk G Wj. Moreover, let wJJ'* 
be the resulting clause of applying the resolution operation on these clauses. In 
this case, if wJJ® is not a tautology, then an edge {uji,ujj) is added to G whose 
weight is defined as: 


w{uJi,UJj) 



(3) 


Notice that in the RES model, clauses are related if the application of the 
resolution operation results in a non-trivial resolvent. Moreover, observe that the 
weight of edges between pairs of clauses is greater when the size of the resolvent 
is smaller. The goal is to make tighter the relations between clauses that produce 
smaller clauses when resolution is applied. 

Consider the following MaxSAT formula where wi : (xi V X 2 ), W 2 : (a ^2 V X 3 ) 
and UJ 3 : {xi V 0 : 3 ) are hard clauses and ui^ : (xi), 025 : (xa) are soft clauses. 
Figures 1(a) and 1(b) illustrate the structure of the graph representation of this 
formula when using the CVIG and RES models. The weights of edges are not 
represented for simplicity. Observe that if the clause luq : (xi V X 2 ) was added to 
the formula, it would not connect to any other clause in the RES graph because 
the only clause containing Xi positively is wi = (xi V X2), but that does not 
connect to ujq due to X 2 appearing negatively and positively in wg ^md wi, re¬ 
spectively. A similar type of analysis is done in blocked clause elimination [13,10] 
— a technique commonly used in formula preprocessing. 

Although resolution-based graphs are not novel [26] and have been used in 
other domains [25], in this paper we propose to enhance the resolution-based 
graph representation by adding weights to edges. Moreover, as far as we know, 
this representation has never been used for MaxSAT solving. 
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4 New Partition-based Algorithm for MaxSAT 


Despite its very good performance in industrial partial MaxSAT instances, the 
MSU3 algorithm (see Algorithm 2) may suffer from two issues: (1) identification 
of unnecessarily large unsatisfiable subformulas and, ( 2 ) a potentially large car¬ 
dinality constraint to be maintained between iterations. In fact these issues are 
related. If an unsatisfiable subformula with an unnecessarily large number of soft 
clauses is encountered early, then an unnecessarily large cardinality constraint 
has to be dealt with through most of the algorithm’s iterations. 

Our approach to tackle these issues is to split the set of soft clauses. The 
goal is that, at each iteration, the algorithm should only consider part of the 
problem, instead of dealing with the whole problem instance in each iteration. 

4.1 Algorithm Description 

Algorithm 3 presents our enhancement of MSU3 with partition of the soft clause 
set. The algorithm starts by partitioning ips into n disjoint sets of soft clauses 
71)72 • ■ •7n (line 1). Observe that several methods can be used to partition Lp^- 
Details of this procedure are discussed later. 

For each set 7 ^, we apply the MSU3 algorithm to the formula iph (lines 2- 
11). As a result, we obtain a lower bound value \i associated with each set of 
soft clauses 7 ^. If the partitioning procedure creates a single partition, then the 
algorithm terminates (line 13). Otherwise, it is necessary to build the solution 
of the MaxSAT instance by merging the different sets of soft clauses. 

The merge process works as follows. At each iteration, two sets of soft clauses 
7 i and 7 ^- are selected to be merged (line 15) and removed from 7 . Let 7 ^ denote 
the union of 7 ^ and 7 ^. Since 7 ^ and 7 ^- are disjoint, we necessarily have that 
Xi -\- Xj is a lower bound for 7 ^. Hence, we can safely initialize Afc = -I- Xj 

(line 17). Next, the lower bound Xk is refined by applying the MSU3 algorithm 
to (pi^ U 7 fc (lines 18-25). When set 7 becomes empty, then all soft clauses were 
merged and the last solution found is an optimal solution (line 27). Otherwise, 
there are still more sets to be merged and 7 ^ is added to 7 (line 29). 

4.2 Partition and Merge of Soft Clauses 

Algorithm 3 can be configured differently depending on two procedures: (1) how 
the set of soft clauses is partitioned (line 1 ) and ( 2 ) how to merge two sets of 
soft clauses (line 15). 

In the partition procedure, our algorithm starts by representing the CNF 
formula as a graph using one of the models described in section 3. Next, we apply 
a community-finding algorithm on the graph representation that maximizes a 
modularity measure [4] in order to obtain a graph partitioning. 

Recently, the use of modularity measures has become widespread when an¬ 
alyzing the structure of graphs, in particular for the identification of communi¬ 
ties [7,23]. In fact, this has already been used in the analysis of SAT instances [1] 
and to improve the initial unsatisfiability-based approach proposed by Fu and 


6 


Algorithm 3: Extended MSU3 Algorithm 


1 

2 

3 

4 

5 

6 

7 
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9 

10 

11 


Input: U 

Output: satisfying assignment to ip 
7 ^ (7i; ■ • ■ ^In) ^ partitionSoft(c^s, pn) 
foreach 7i ^ 7 do 
(VS,AO^(0,O) 

(st, (fie, ^ SAT((^h U 7i) 
while St = UNSAT do 

foreach ut € {^pc H pa) do 

Vfl ■(— Vr U {r} // r is a new variable 

ujR ■<— OJ U {r} // uj was not previously relaxed 

L 7i ^ (7i \ U {u)b} 

\i i — Ai + 1 

[ (st, pcv)^ SAT((^;, U 7 i U {CNF(E,gyj r < Ai)}) 


12 if I 7 I = 1 then 

13 ^ return!/ // no partitions were identified 


14 while true do 

15 (7i, 7j)-f—selectPartitions(7) 

16 7^7\{7i,7j} 

17 (7,, 1/|, Afc) ^ (7! U 7 i, VS U VS, Ai + A,-) 

18 (st, pc, v) <r- SAT(vJh u 7fc U {CNF(5S^gv'= ’’ ^ Afe)}) 


19 

20 
21 
22 

23 

24 

25 


while St = UNSAT do 

foreach u € (pc H ps) do 

// r is a new variable 
cjj? ■<—CJ U {r} // uj was not previously relaxed 

|_ 7fc ^ (7fc \ {^^}) U {tufl} 

Afc i — \}z “h 1 

(st, pc, y) ^ SkT{ph u 7fe U {CW(Y^^^yk r < Xk)}) 


26 

27 

28 
29 


if 7 = 0 then 
return u 
else 

|_ 7^7U{7fc} 


Malik [6,18]. The purpose of the modularity measure is to evaluate the quality of 
the partitions, where vertices inside a partition should be densely connected and 
vertices assigned to different partitions should be loosely connected. However, 
finding a set of partitions with an optimal modularity value is computationally 
hard [5]. In our implementation, we use the approximation algorithm proposed 
by Blondel et al. [4]. 

At each iteration in Algorithm 3, two partitions are selected to be merged. 
One can devise several different criteria to select and merge the partitions of 
soft clauses. In early attempts, the merge process was sequential [18]. Given n 
partitions 71,72 ... 7 „, at iteration i (i < n) of the algorithm, the first i partitions 
7 i 172 • ■ • 7 i were merged sequentially. 


7 












(a) Sequential merging (b) Balanced merging 


Fig. 2. Examples of merge processes 


Figure 2(a) illustrates the sequential merging procedure. Observe that the 
sequential merging process is not balanced. This results in an early growth of 
the identified subformulas and, as a result, an early growth of the cardinality 
constraints to be maintained at each iteration of the algorithm. 

In this paper, we propose a weighted balanced merge procedure that depends 
on the strength of the graph connections between partitions. The goal is to delay 
having to deal with a large number of soft clauses, until the latter iterations of 
the algorithm. Figure 2(b) illustrates the weighted balanced merging procedure. 

Let G = {V, E) denote an undirected weighted graph where V is the set of 
vertices and E the set of edges. Let w : if —>■ K be a weight function for each edge 
in the graph. The community-finding algorithm identifies a set of communities 
C = {Ci,C 2 , ■ ■ ■ ,Cn} where every vertex u G V is assigned to one and only one 
community in C. Hence, since in both CVIG and RES model there is a node 
for each propositional clause, one can build the partitions in a straightforward 
manner. For each community Ci with vertices representing soft clauses, there is 
a partition 7 ^ containing the respective soft clauses. 

Based on the graph representation, one can define the strength of the connec¬ 
tion between partitions. Let dij denote the strength between partition 7 ^ and 7 ^. 
One can define dij based on the weight between the vertices of their respective 
communities Ci and Cj in the graph. Hence, dij can be defined as follows: 

w{u,v) (4) 

uGCiAvGCj 

Considering that the graph is undirected, we necessarily have that dij = dji. 

Given an initial set 7 of n partitions 71,72 • ■ •, 7 n > our algorithm applies a 
greedy procedure that pairs all partitions 7 i and 7 ^ from 7 to be merged, starting 
with the pair with largest dij. After pairing all partitions in the initial set, we 
perform the same procedure to the next n /2 partitions that result from the 
initial merging iterations. This is iteratively applied until we only have a single 
partition (see Figure 2(b)). 

Observe that if partitions 7 ^ and 7 ^- are merged into a new partition 7 ^, 
then the connectivity strength dki between jk to another partition 7 ; is given by 





































dki = dii + dji- This follows from the fact that the communities in the graph are 
disjoint. 

Finally, we would like to reference other solvers that split the set of soft 
clauses by identifying disjoint unsatisfiable subformulas [8,21]. However, there are 
major differences with regard to our proposed approach. First, our solver takes 
advantage of an explicit formula representation to split the set of soft clauses, 
instead of using the unsatisfiable subformulas provided by the SAT solver. More¬ 
over, in our solver, the merge process is also guided by the explicit representation 
of the formula. 

Furthermore, in solvers where disjoint unsatisfiable subformulas are iden¬ 
tified [ 8 , 21 ], the split occurs on the cardinality constraints at each iteration. 
However, each SAT call still has to deal with the whole formula at each itera¬ 
tion. In Algorithm 3, the SAT solver does not have to deal with all soft clauses 
at each iteration, but only after the final merge step. 

4.3 Algorithm Analysis 

In this section a proof sketch of the correctness Algorithm 3, as well as an analysis 
on the number of SAT calls is presented. 

Proof. As mentioned in section 2 , we assume the set of hard clauses (ph is satis- 
fiable. Otherwise, the MaxSAT formula is unsatisfiable. This can be verified by 
a single SAT call on ph before applying Algorithm 3. 

For the proof we adopt the following notation. For some set processed in 
Algorithm 3, we write 7 /^ C for the set of clauses that were relaxed in the 
algorithm (but clauses in 7 /^ do not contain the relaxation variables). We will 
prove by induction the invariant that ph U 7 /^ cannot be satisfied unless at least 
Xi clauses are removed from 7 /^. The induction hypothesis is satisfied trivially 
at the beginning of the algorithm as each Xi is initialized to 0 . 

Consider the case where Xi is augmented by 1 when Ph^li^{J2reV' ’’’ — 
is unsatisfiable. Let pc be the obtained unsatisfiable subformula from the SAT 
call, let Pq C p^ be the soft clauses of pc that appear as relaxed in 7 ^ and let 
Pq = Ps C] Pc be the rest of the soft clauses in the unsatisfiable subformula 
(not yet relaxed). From induction hypothesis ph U p^ cannot be satisfied unless 
at least Xi clauses are removed from p^ C 7 ^^. Since pc is an unsatisfiable 
subformula, it is impossible to satisfy ph U p^ U Pq by removing Xi clauses 
from p^. Now we need to also show that it is impossible to satisfy ph U 7 /^ U Pq 
by removing Xi clauses from 7 /^ U Pq (this is the new set of relaxed clauses). 

Let us assume for contradiction that it is possible to satisfy 7 /^ U p^ by 
removing some set of clauses ^ s.t. j^j = Xi. To show the contradiction we 
consider two cases: (1) f Q ji and (2) ^ ^ 7 .^. Case (1) yields an immediate 
contradiction as we would have not obtained unsatisfiability in the SAT call 
as it would be possible to satisfy ph U 7 /^ by removing Xi clauses from 7 /^. For 
case ( 2 ) consider that there is a clause ui G ^ s.t. oj is not yet relaxed, i.e. uj ^ 7 ^. 
This means that ph U 7 /^ is satisfiable after removing less than Xi clauses, which 
is a contradiction with the induction hypothesis. 
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To show that the invariant is preserved by the merge operation, we observe 
that any merged 7 ^ and 7 ^ are disjoint and therefore so are 7 /^ and 7 ^. In order 
to satisfy both (^/jUy/^, must be satisfied. Consequently, 

at least Xi + Xj clauses must be removed from (yf U yj^). 

□ 

Finally, we note that the number of SAT calls performed by Algorithm 3 
is larger than the MSU3 algorithm. Observe that the number of unsatisfiable 
SAT calls is the same for both algorithms. Let A be the number of unsatisfiable 
soft clauses at any optimal solution of the MaxSAT instance. In this case, both 
algorithms perform A unsatisfiable SAT calls. However, while MSU3 performs 
only one satisfiable SAT call. Algorithm 3 performs 2n —1, where n is the number 
of identified partitions (line 1 ). 


5 Experimental Results 

In this section we compare different configurations of Algorithm 3 with the top 3 
non-portfolio solvers of the MaxSAT 2014 Evaluation’s partial MaxSAT category. 
The top 3 were Open-WBO’s MSU3 incremental algorithm [17,16], EvaSOOa [22] 
and MSCG [20] . The new partition-based algorithm is also implemented using 
the Open-WBO framework^. 

The algorithms were evaluated running on the power set of the partial MaxSAT 
industrial instances of the MaxSAT evaluations of 2012, 2013 and 2014. For each 
instance, algorithms were executed with a timeout of 1800 seconds and a memory 
limit of 4 GB. Similar resource limitations were used during the last MaxSAT 
Evaluation of 2014. These tests were conducted on a machine with 4 AMD 
Opteron 6376 (2.3 GHz) and 128 GB of RAM, running Debian Jessie. 

Table 1 presents the number of instances solved by each algorithm, per in¬ 
stance set. Besides MSU3, EvaSOOa and MSCG, results for the best 4 config¬ 
urations of the partition-based enhanced MSU3 algorithm are shown. S-CVIG 
applies the sequential merging of partitions using the CVIG graph model. S- 
RES also applies sequential merging, but using the RES graph model. W-CVIG 
and W-RES apply the weighted balanced merging of partitions, using the CVIG 
and RES graph models, respectively. Note that all our implementations are fully 
incremental, i.e. only one instance of the SAT solver is created throughout the 
execution of the proposed algorithm. As with the MSU3 implementation on 
Open-WBO, we take advantage of assumptions usage at each SAT call and in¬ 
cremental encoding of cardinality constraints [16] . 

Results from Table 1 show that all variants of the partition-based algorithm 
are competitive with the remaining state of the art algorithms. However, overall 
results clearly show that W-RES outperforms all remaining algorithms, since it 
is able to solve more instances in total. Moreover, results for the configurations 
of partition-based algorithm also show that weight-based balanced merging of 
partitions is preferable to sequential partitioning. 

^ Available at http://sat.inesc-id.pt/open-wbo/ 
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Instance Group 

Total 

MSU3 

EvaSOOa 

MSCG 

S-CVIG 

S-RES 

W-CVIG 

W-RES 

aes 

7 

1 

1 

1 

1 

1 

1 

1 

atcoss/mesat 

18 

11 

11 

4 

11 

1 

11 

11 

atcoss/sugar 

19 

12 

11 

4 

12 

3 

12 

12 

bcp/fir 

59 

59 

55 

59 

56 

44 

51 

51 

bcp/hipp-yRal/simp 

17 

16 

16 

16 

16 

16 

16 

16 

bcp/hipp-yRal/su 

38 

35 

34 

33 

34 

34 

35 

33 

bcp/msp 

64 

26 

37 

29 

23 

41 

27 

42 

bcp/mtg 

40 

40 

40 

40 

40 

40 

40 

40 

bcp/syn 

74 

43 

48 

47 

47 

48 

46 

49 

circuit-trace-compaction 

4 

4 

4 

4 

4 

3 

4 

4 

close-solutions 

50 

48 

48 

46 

40 

32 

40 

45 

des 

50 

42 

41 

41 

49 

48 

50 

48 

haplotype-assembly 

6 

5 

5 

5 

5 

5 

5 

5 

hs-tinietabling 

2 

1 

1 

0 

1 

1 

1 

1 

mbd 

46 

45 

42 

43 

44 

45 

45 

45 

packup-pms 

40 

40 

40 

40 

40 

40 

40 

40 

pbo / mqc / nencdr 

84 

84 

84 

84 

84 

84 

84 

84 

pbo / mqc / nlogencdr 

84 

84 

84 

84 

84 

84 

84 

84 

pbo/routing 

15 

15 

15 

15 

14 

15 

15 

15 

protein ins 

12 

12 

8 

12 

12 

12 

12 

12 

tpr/Multiple path 

48 

48 

44 

42 

48 

48 

48 

48 

tpr/One path 

50 

50 

50 

50 

50 

50 

50 

50 

Total 

827 

721 

719 

699 

715 

695 

717 

736 


Table 1 . Experimental evaluation of Open-WBO’s MSU3 algorithm, EvaSOOa, MSCG 
and 4 different configurations of the partition-based algorithm. 


Considering that MSU3 is our base solver, most gains occur in instance sets 
bcp/msp, bcp/syn and des. While in the bcp/syn and des instance sets, all 
partition-based configurations perform better, in bcp/msp the resolution-based 
graph partitioning allowed a significant performance boost. 

Figures 3(a) and 3(b) compare the results of S-RES and W-RES on the des 
and bcp/msp instance sets. In the des instances, the run time of sequential merg¬ 
ing is slightly better, despite solving the same number of instances. Nevertheless, 
in the bcp/msp instance set the weight-based balanced merging used in W-RES 
clearly outperforms the sequential merging approach used in S-RES. 

In Figures 4(a) and 4(b) we compare MSU3 and W-RES on the same bench¬ 
mark sets. It can be observed that W-RES performs much better in these in¬ 
stances. In the des instance set, there are some instances where W-RES is not as 
fast, since there is some time spent in finding partitions and additional SAT calls. 
We note that there is always some time spent in building the graph, applying 
the community finding algorithm and splitting the set of soft clauses. However, 
this partitioning step is usually not very time consuming. Nevertheless, W-RES 
is able to scale better and solve more instances. In the bcp/msp instances, the 
proposed techniques allow W-RES to be much better than MSU3, as well as all 
other algorithms tested. 

Resolution-based graph models performed worst in the bcp/fir category. 
It was observed that the overall modularity values obtained for the resolution- 
based graphs were low in this particular instance set. As a result, the partitioning 
obtained for S-RES and W-RES in bcp/fir instances is not as meaningful as for 
other instance sets. Hence, when this occurs, it can deteriorate the algorithm’s 
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(a) des instance set 



(b) bcp/msp instance set 


Fig. 3. Comparison between run times of S-RES and W-RES on des and bcp/msp 
instance sets 


performance, since the partition-based algorithm performs more SAT calls than 
MSU3. 

Finally, Figure 5 shows a cactus plot with the run times of all algorithms 
considered in the experimental evaluation. Here we can observe that S-RES is 
much slower than W-RES, clearly showing the effectiveness of the newly pro¬ 
posed weight-based merging. Overall, W-RES clearly outperforms the remaining 
algorithms, being able to solve 700 instances in 300 seconds or less. 


6 Conclusions and Future Work 

In this paper we exploit resolution-based graph representations of CNF formulas 
in order to develop a new state of the art algorithm for MaxSAT. In the proposed 
approach, soft clauses are initially partitioned in disjoint sets by analyzing the 
formula structure. The partitioning process is attained by applying a community¬ 
finding algorithm on weighted resolution-based graphs. Next, at each iteration 
of the algorithm, partitions are merged using structural information from the 
graph representation until an optimal solution is found. 

The proposed approach is novel in many aspects. First, the use of a resolution- 
based graph representation allows to better model the interaction between clauses. 
Furthermore, instead of applying a sequential merging process, the graph repre¬ 
sentation is also used in a weight-based balanced merging procedure. Moreover, 
since the algorithm does not have to deal with the whole formula at each iter¬ 
ation, smaller unsatisfiable cores are identified. As a result from this process, 
smaller cardinality constraints are encoded into CNF at each iteration, thus 
improving the algorithm’s performance. 
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(a) des instance set 


(b) bcp/msp instance set 


Fig. 4. Comparison between run times of MSU3 and W-RES on des and bcp/msp 
instance sets 



Fig. 5 . Cactus plot with the run times for MSU3, EvaSOOa, MSCG, S-CVIG, S-RES, 
W-CVIG and W-RES. 


Experimental results obtained in industrial partial MaxSAT instances clearly 
show the effectiveness of the proposed algorithm. As a result, our solver improves 
upon the best non-portfolio solver from the 2014 MaxSAT solver evaluation. 
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The source code of the new solver will become available as part of the Open- 
WBO framework. This will allow the research community to build upon the 
current work to further improve MaxSAT solving. 

As future work, we propose to extend the proposed approach for weighted 
MaxSAT solving. Moreover, different model representations of CNF formulas are 
to be tested, as well as new techniques for building and merging partitions of 
soft clauses in MaxSAT formulas. Furthermore, the proposed techniques are not 
exclusive to MSU3 and can also be integrated into other MaxSAT algorithms. 
Additionally, these techniques can also be applied to other extensions of SAT. 
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